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AN EXPERIMENTAL/ANALYTICAL PROGRAM TO ASSESS THE UTILITY OF 
LIDAR FOR POLLUTION MONITORING 


By 

Frank S. Mills 1 , Robert J. Allen 2 , and Carolyn F. Butler 2 


1 . INTRODUCTION 

The aim of the research performed under this grant was the development 
and demonstration of lidar techniques for the remote measurement of atmos- 
pheric constituents and transport processes in the lower troposphere. 
Particular emphasis was given to techniques for monitoring S0 2 and partic- 
ulates, the principal pollutants in power plant and industrial plumes. 

Data from a plume dispersion study conducted in Maryland during September 
and October 1976 were reduced, and a data base was assembled which is 
available to the scientific community for plume model verification. (The 
mobile ruby lidar system used in the study was constructed under a previous 
grant, reference 1). A UV Differential Absorption Lidar (DIAL) was built, 
and preliminary testing was done. Work on the development of an IR DIAL 
system proceeded. 

This report will describe in detail the accomplishments in all areas 
mentioned above. Section 2 describes a study of the effect of transient 
digitizer errors on DIAL measurements conducted by Frank S. Mills of Old 
Dominion University and Edward V. Browell of NASA Langley Research Center. 

A paper describing this study was presented at the Eighth International 
Laser Radar Conference in June 1977 (ref. 2). This study is pertinent to 
all DIAL systems including the UV and IR DIAL systems. 

Section 3 contains a description of the UV DIAL system. The system 
was designed and constructed by Robert J. Allen and Lawrence Perillo of 
Old Dominion University working with NASA LaRC personnel using a double- 

1 Research Assistant Professor of Physics and Geophysical Sciences, Old 
Dominion University, Norfolk, Virginia 23508. 

2 Graduate Research Assistant, Department of Physics and Geophysical 
Sciences, Old Dominion University, Norfolk, Virginia 23508. 


pul so laser concept developed by Dr. Browell at NASA LaRC. Dr. Mills was 
responsible for data acquisition system hardware and software for the UV 
DIAL. A paper describing the UV DIAL system was presented at the Eighth 
International Laser Radar Conference in June 197? (ref. 3). 

Section 4 briefly describes progress made on the IR DIAL system, 

Progress was limited since the laser delivery was delayed due to problems 
encountered by the laser manufacturer. Much of the data acquisition 
system software developed for the UV DIAL system will also be applicable 
to the IR DIAL system. 

Section 5 is a report of the plume dispersion study which was conducted 
during September and October 1976 including a description of the mobile 
lidar system, the experimental design, the data reduction process, and the 
resulting data base. This study was sponsored on a cooperative basis by 
NASA and the State of Maryland Power Plant Siting Program with participa- 
tion by a number of other organisations including Old Dominion University, 
Wyle Laboratories, Martin Marietta Research Laboratories, Environmental 
Measurements Incorporated, Engineering-Test Services, and the Potomac 
Electric Power Company. Old Dominion University personnel involved in the 
study were Dr. Mills and Carolyn F. Butler, A paper describing the study 
was presented at the Eighth International Laser Radar Conference in 
June 1977 (ref. 4). 

During the grant period, Ms. Butler completed a two-dimensional modeling 
study of the stratospheric transport of volcanic aerosols. This work has 
been reported elsewhere (refs, 5, 6), 

Also during the grant period, Neale Mayo, a graduate research assistant 
in the Department of Physics and Geophysical Sciences at Old Dominion Univer- 
sity, worked on the development of lidar systems at NASA LaRC. Much of his 
work is reported here. 

2. THE EFFECT OF TRANSIENT DIGITIZER ERRORS ON DIAL MEASUREMENTS 

This study was prompted by errors encountered in analyzing atmospheric 
returns from the near-IR DIAL system at NASA Langley Research Center which 
is used to measure vertical water vapor profiles. Figure 1 shows an example 
of the problem encountered. An exponentially decaying test signal was fed 



through both channels of the data system, and a concentration profile was 
calculated in exactly the same way as for atmospheric returns. The large 
excursions at 1.0 and 2.0 km are caused by amplifier switching transients. 
The anomalies at 0.35 and 1.5 km are the ones which motivated this study. 

In order to understand these anomalies it is necessary to describe how 
the concentration is calculated, and then to describe how the transient 
digitizer operates. 

The DIAL concept is described in detail in section 3. Therefore only 
those aspects pertaining to this study will be discussed here. 

The average concentration of the gas of interest in the range between 
ranges and R 2 is given by 
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where a. (X ) - a . (X is the difference in the absorption cross section 
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of the gas of interest at the on and off wavelengths; P„ _ fR.) is the 
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power received at the wavelength from range R^» and P r,off tr V is the 
power received at the off wavelength from range R^. Before the concen- 
tration is calculated, the received signals are usually smoothed by 
performing a running average. The averaging interval is typically the 
same as the range cell length R 2 - although a smaller interval could 
just as easily be used. 


One of the transient digitizers studied was the Biomation model 8100. 
This unit is capable of sampling an analog input at sample intervals as 
short as 10 ns, converting the analog samples to digital form, and storing 
the samples in a 2048-word memory. The digital resolution of the unit is 
8 bits or 1 part in 256. For sample intervals as short as 100 ns the unit 
has the additional capability of accepting 2 independent analog inputs, 
which are designated channel A and channel B. In this mode of operation 
the analog inputs are sampled alternately, and the wave forms are stored in 
memory so that the even-numbered points contain the digital equivalent of 
the channel A input and the odd-numbered points contain the digital equiv- 
alent of the channel B input. For each channel the input signal is fed 



through an input attenuator and a fixed gain video amplifier. The output 
of the video amplifier is then fed through a folding amplifier or full 
wave rectifying amplifier and then to the analog to digital converter. 

The folding amplifier passes positive signals unchanged and inverts nega- 
tive signals. A sign bit is generated if the input was initially negative 
and inverted by the folding amplifier. Thus the folding amplifier produces 
the most significant or "$ign M bit of the analog to digital conversion and 
reduces the dynamic range requirement for the remainder of the conversion 
by a factor of two. The remaining seven bits of analog to digital conversion 
are produced by a parallel comparator type of converter. 

The primary difficulty encountered with this type of transient digitizer 
is a discontinuity in the digitized wave form at the point where the input 
signal crosses zero. This is the source of the anomalies shown in figure 
1. There is an adjustment in the digitizer, the crossover adjustment, 
which can be set to eliminate the discontinuity, but the adjustment is 
very delicate and temperature sensitive, particularly when the unit is 
operating in the two-channel mode. Three different Biomation 8100 digitizers 
were investigated, and they all exhibited the same behavior, which indicates 
that the problem is inherent in the design and not just a peculiarity asso- 
ciated with one unit. 

Figure 2 shows schematically how the calculated concentration is affected 
by the two different types of crossover misadjustment . It should be pointed 
out that if the discontinuity is of the same kind and magnitude and occurs at 
the same point in time for both the on and off wavelength return, there will 
be no anomaly in the calculated concentration, since the ratio of received 
powers in the argument of the logarithm in equation (1) remains unchanged. 

A number of tests were performed on the Biomation 8100 using a decaying 
exponential signal ta the input to simulate lidar returns. In the first 
series of tests, the same signal was fed to both channel A and channel B. 

The crossover adjustment for each channel was carefully adjusted to eliminate 
the discontinuity at the zero crossing point. This was verified by observing 
the digitized wave form for each channel as reproduced by the computer. 

One hundred wave forms on each channel were averaged and the concentration 
calculated assuming that a water vapor measurement was being made. For 
this test the calculated concentration showed no anomaly at the zero 
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crossing point. The second test used the sane data except that the channel 
A data was shifted with respect to channel B by an amount which would be 
equivalent to S90 m before the concentration calculation was made. In this 
case, since the zero crossings occur at different times, the anomaly asso- 
ciated with the crossover misadjuatment appeared even though every effort 
was made to adjust the crossover properly. 

It was observed while making the crossover adjustment that the adjust- 
aent appeared to be most unstable and sensitive when the positive and nega- 
tive parts of the wave form were correctly matched. It was therefore decided 
to adjust the crossover for each channel so that there was a small gap in the 
digitized wave form similar to the left-hand side of figure 2 and correct for 
the aisadjustment analytically. First a slowly decaying signal was fed to 
each channel, and the average size of the gap was measured for each channel 
for 100 input wave forms. Next the first test was repeated, and the gap 
previously measured for each channel was subtracted from the portion of each 
input wave form which was greater than zero. As before, 100 wave forms on 
each channel were averaged and the concentration calculated with channel A 
shifted 390 m with respect to channel B. The result was an anomaly of the 
type shown in the right-hand side of figure 2, implying that the gap correc- 
tion was too large. 

The same input signals were reprocessed using a slightly smaller gap 
correction for channel A. The result is shown in figure 3. The feature 
between the third and fourth tick marks is due to the channel A midscale 
problem, and the feature between the fourth and sixth tick marks is due to 
the channel B midscale problem. Examining the channel A feature more care- 
fully, it can be seen that there is gap-induced anomaly (as shown in the 
left side of figure 2) which arises because the gap correction is now too 
small, superimposed on a broader anomaly, caused by a flat spot in the 
return (as shown in the right side of figure 2). Since the crossover 
adjustment has been set to produce a gap in the return, it cannot also 
produce a flat spot in the return. Therefore, the flat spot is caused by 
something else. The most likely source is a small nonlinearity in the 
folding amplifier when the input signal is near zero. 

The tests just described have demonstrated problems which are asso- 
ciated with the point where the input signal crosses zero. The first is 
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the sensitivity of the crossover adjustment, The second is a slight non- 
linearity* probably in the folding amplifier. The crossover adjustment 
problem, while annoying, is not insurmountable, since it can be corrected 
analytically. 

The second problem, however, is more serious and would be quite 
difficult to correct* It should be pointed out that the nonlinearity at 
the zero crossing point is a very small effect and in most applications is 
not detectable. However, the differential absorption calculation is very 
sensitive to nonlinearities. Therefore, for differential absorption appli- 
cations, the input signals should be restricted so that they are always 
positive or always negative, thus avoiding the pioblems at the zero crossing 
point. This reduces the dynamic range from eight to seven bits so it may 
be necessary to average more returns to get acceptable accuracy. 

The tests described above were repeated with the input signals strictly 
greater than zero. The calculated concentration showed no systematic 
anomalies. 

The second transient digitizer studied was the Biomation model 1010. 
This digitizer can sample an analog input wave form at intervals as short as 
100 ns, convert the analog samples to digital form, and store the samples 
in a 2050-word memory. The digital resolution of this unit is 10 bits or 
1 part in 1024. The input signal passes through the input amplifier which 
consists of two identical stages, each preceded by an attenuator section 
to provide the various input ranges, and then to the analog to digital 
converter. The analog to digital converter is a dual -rank type, or two 
parallel comparator converters in series. The first parallel comparator 
converter determines the five most significant bits of the digital result, 
and the second parallel comparator converter provides the five least 
significant bits. 

The performance of the Biomation 1010 digitizer was investigated using 
a decaying exponential signal to simulate a lidar return. The model 1010 is 
a single channel device, so two-channel operation was simulated by treating 
the even-numbered words of the digitized signal as the on-line return and 
the odd-numbered words as the off-line return. As before, 100 wave forms 
were averaged and the concentration calculated with the on-line data shifted 
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with respect to the off-line data so that any systematic errors would show 
up. The result showed no systematic error which could be attributed to the 
transient digitiser. 

There are two primary conclusions to be drawn from this study. The 
first is that for sampling intervals of 100 ns or greater the Biomation 
model 1010 or a similar unit should be used for differential absorption 
measurements. Second, if sampling intervals shorter than 100 ns are 
required, the Biomation model 3100 can be used, but input signals must be 
either strictly positive or strictly negative in order to avoid the error 
at the zero crossing point. 


3. UV DIAL SYSTEM 

The UV DIAL system is being developed for the remote investigations of 
SO2, O3, and aerosols present in power plant and urban plumes. The data 
obtained will be used together with that from other sensors in dispersion, 
transformation, and deposition studies. It is important to know what 
happens to the SO2 because of the potential health and environmental 
effects of sulfur compounds. In addition, these studies will allow the 
reassessment of the accuracy of current plume models in predicting the 
temporal and spatial changes in pollutant concentrations and compositions. 
It is important to study O3 because of its role in photochemistry and its 
potential adverse effect: cn health and vegetation in large concentrations 
(greater than SO- 100 ppb) . 

Development efforts during the past year have been primarily with the 
laboratory SO2 1 dar system operating near 300 nm. Future plans are to 
"harden” the system for use in the field and aboard an aircraft, and to 
develop the system capabilities for O3 investigations around 283 nm. 

The system uses a double-pulse frequency-doubled Neodymium-doped 
yttrium-aluminum garnet (Nd:YAG) laser as a pump for two independent dye 
lasers. This concept, developed by Dr. Browell at NASA/LaRC (ref. ?), 
allows the transmission of two closely spaced sequential laser pulses 
at two different wavelengths (300.2 nm and 299.5 nm) for SO2 differen- 
tial absorption observations. In addition, the system produces an output 
at 500.4 nm for aerosol backscattering observations. The DIAL concept 



will new be described qualitatively. This section will be completed with a 
description of the IIV DIAL system 

3.1, UV DIAL Concept 

(See reference 1.) The Differential Absorption Lidar (DIAL) has become 
a primary method for remote measurement of gases in the troposphere. To 
accomplish this, pulsed laser radiation at two wavelengths is transmitted 
into the atmosphere. The two wavelengths are selected so that one, called 
the on wavelength, is absorbed by the gas of interest, while the other, 
called the off wavelength, is not. The backscatter return signal as a 
function of range at each wavelength is collected by an optical receiver. 

Now define a range cell as the distance from Rj to R 2 . The average concen- 
tration of the absorbing gas in the range cell may be determined from the 
on and off wavelength returns at Ri and R 2 . The situation is analagous to 
a dual-beam spectroscopy experiment where the off wavelength return corre- 
sponds to the reference beam, the on wavelength corresponds to the sample 
beam, and the range cell corresponds to the sample cell. The transmittance 
of the absorbing gas in the range cell is then just the ratio of the on 
wavelength returns at R 2 and Ri divided by the ratio of the off wavelength 
returns at R 2 and Rj. 

The quantitative expression for absorber concentration can be found 
using the lidar equation: 

P r (R,X,t) - 1C i eCR.X.t) 

• exp |-2/ R £i jc (r,X,t) ♦ N A (r,t)o A (X,P,T) 

* N int (r ’ t >'’int CX - P - T ’] dr ) 

Here P is the received power as a function of range R, wavelength X , and 
r 

time t; K is a constant; L is the laser pulse length; A is the receiver area; 

P (X), the transmitted power, is a function of wavelength; S is the back- 
o 

scatter cross section as a function of range, wavelength, and time; is 
the extinction coefficient for scattering as a function of range, wavelength 
and time; is the concentration of the absorbing gas of interest as a 
function of range and time; is the absorption cross section of the gas of 
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interest as a function of wavelength, pressure, P, and temperature, T; N^ nt 
is concentration of any interfering absorbing gas; and is the absorp- 
tion cross section of the interfering gas. 

A number of assumptions are made to simplify t' ■* analysis. The scat- 
tering terms are assumed equal for the on wavelength and the off wavelength 
returns, that is 8(R,Aoff,t off) ■ B(R,Aon,t on) and at (R,Ao£f,t off) * 

a (R.Aon.t on). Also, the * jncentrations of the absorbing and interfering 
gases are assumed to remain constant between the off wavelength and on wave- 
length returns, that is, N A (R,t off) * N A (R,t on) and N int ( R » t * N int 
(R,t on). The on and off wavelength absorption cross section for the gas 
of interest o A (Aoff,P,T) and J A (Aon,P,T), are assumed to be known and 
constant between Rj and R 2 . If Oj (Xoff,P,T) is not equal to o int 
(Aon,P,T), then both the concentration and absorption cross section of the 
interfering gas must be known or must be determined by a separate experiment. 

With the above assumptions the average concentration of >L of the gas 
of interest in the range cell between R] and R 2 is 
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The small separation in time (less than 100 us) of the on and off wave- 
length pulses in the UV DIAL transmitter minimises errors in the measured 
concentration which might arise from time variations of 8 and a s{; . The 
time separation also allows the receiver electronics to be time multiplexed, 
with the same photomultiplier and transient digitiser being used for both 
the on and off wavelength returns, thus eliminating another potential source 
for systematic error. 

The DIAL technique was first used by Schotland (ref. 8) to measure 
water vapor. Subsequent remote measurements of NO- (ref. 9) , S0 2 and O3 
(ref. 10) have been made by other investigators using this method. 


3.2. Description of the UV DIAL System 

The UV DIAL System consists of a trilaser transmitting system, a dual- 
photomultiplier multiplexed receiver, data processing subsystem, and 



associated timing/control logic. The transmitter uses an International 
Laser System (ILS) double-pulsed frequency-doubled Nd:YAG laser as a pump 
for two dye lasers, each tuned to a slightly different wavelength. The 
ILS laser, designed to NASA specifications, produces a pulse-pair (pulse A 
and pulse S) transmitted at a wavelength of 532 nm sequentially at time 
t and (t Q + At) where 90 At 99 ps. The energy contained in each pulse 
is about ICO aJ with a maximum prf of S Hz. 

Pulses A and B are directed either to an off-line or an on-line dye 
laser by an optical switch consisting of a Pockels cell, dielectric 
polarizer and half-wave plate (see figure 4). If the applied voltage to 
the Pockels cell crystal is zero, the horizontally polarized pulse from 
the ILS laser will pass through the cell unaffected and be directed to 
the off-line dye laser by the dielectric polarizer. When a high-voltage 
pulse is applied to the Pockels cell, the transmitted polarized light is 
rotated 90®, thereby being directed by the dielectric polarizer to the 
on-line cavity. A half-wave plate is inserted to rerotate the light back 
to its original horizontal polarization. A switch S-l has ueen added to 
the Pockels cell control electronics that will allow the option to open 
the Pockels cell at time t or at time (t + At), thereby routing pulse 
A either to the on-line or the off-line dye laser respectively. This is 
summarized in table 1. 


Table 1. Optical switch pulse routing, off- or on-line dye laser. 


Switch S-l 

Pockels Cell 


Pulse 

Directed to 

Position 

Open at Time: 

On-line 

Laser 

Off-line Laser 

A on-line 

t 

0 

A 


B 

B on-line 

t + At 

B 


A 

0 





The components in both dye lasers are identical. These consist of a 
beamsplitter, oscillator and amplifier dye cells, an oscillator cavity 
output mirror, and a 79-grove/mm echelle grating. The beamsplitter divides 
the green wavelength energy transmitting about 45 mJ to pump the oscillator 
cell and reflecting 105 mJ to pump the amplifier cell. The dye is 
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circulated through all four cells in series* and consists of a 1.5 * io' 4 
molar concentration of Kiton Red S dye dissolved in water with 3 percent 
Ammonyx 10, which makes the solution slightly basic. The dye cavities 
are presently tunable over a range from about S85 nm to 615 run by the 
echelle gratings. The on-line cavity is nominally tuned to 600.4 nm and 
the off-line to 599 nm for SOg observations. The dye will be changed and 
probably circulated through only two cells in the on-line laser for optimum 
operation at 566 nm for the future 0 3 observations. By using the dye 
amplifiers, the oscillator output energy is increased to 25 mj with a line 
width of 0.04 nm. 

The horizontally polarized output from each of the two dye lasers is 
individually directed through two temperature- tuned ADA frequency doubling 
systems. The ADA crystals double the dye laser wavelength into the UV 
with about 10 percent conversion efficiency producing about 2.S mj of UV 
3 ,?*-' about 22,5 mj of the fundamental. A dichroic beamsplitter following 
tto doubler passes the fundamental for propagation into the atmosphere and 
reflect* the second harmonic through the beamsplitter/mirror arrangement 
shown in figure 4. 

The beamsplitters in the UV primary paths direct a small portion of 
the UV energy to a SO 2 absorption cell assembly. Two adjustable irises 
are mounted on the two input apertures (for pulses A and B) located on the 
side and end of a black cover box. A beamsplitter, absorption cell, 
diffusers, UV neutral density filter, and two photodiodes are mounted 
inside the area covered by the box. The arrangement of these components 
is shown in figures 5 and 6 The UV energy from each pulse is directed by 
the beamsplitter so that approximately one-half the energy is directed 
through the absorption ceil to the photodiode labeled g^, while the 
remaining energy bypasses the cell and is directed to the photodiode 
labeled g„ B . This arrangement together with the calibration and normal!- 

Up 

* Studies were conducted to determine the optimum dye concentration as a 
function of power out of the dye oscillator cavity and amplifier. It 
was found that there was only a slight decrease (6%) in output power by 
using the same dye concentration and circulating the dve in series 
within all tcur dye culls. It was felt that the simplified dye-circu- 
lating system v. ..‘ranted the sacrifice in power at this time. 
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lation process discussed below provides s means for measuring the absorp- 
tion coefficient, a(X), with each lidar shot and accounting for uncer- 
tainities in wavelength, X. 


3.3. Absorption Cell Theory 

Figure S shows the optical path when the UV off-line pulse enters the 
side of the absorption cell assembly; Figure 6 shows the on-line pulse 
entering the end of the absorption cell assembly. These two could be inter- 
changed by the t - (t ♦ t) switch on the Pockels cell control unit as 
discussed earlier in this report. 

From Figure 5 we can write, 


V 3 (t) ■ I off Ct) R off T c T (X) T R S DT G AT * dc offset 


(•*) 


and 
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Likewise Figure 6 produces 
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where I(.t) and I ( t * ) are the intensities of the UV laser beams, after 
passing the iris, R represents the reflectance and T the transmittance of 
the beamsplitter, and T,(\), T g (X), and T^CX) represent the transmittance 
of the empty cell, absorbing gas, and filters, respectively. The terms 
gp and are the transfer characteristics KG(jw) of the detectors 
and amplifiers (quad- integrator} respectively. The dc term is the 

combined dc offset from the battery and drift within the amplifier, which 
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is minimised by capacitance coupling and calibration of the quad -integrator. 
The second subscripts T and B represent the path through and the path by- 
passing the cell respectively. The transmittance of the gas within the cell 
can be further defined by 


T (X) • e 
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•oLP 
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where L is the path length of the absorbing gas in cm, P is the gas pressure 
in atmospheres, and o is the absorption coefficient in (atra-cm)'" 1 . 

In order to obtain an expression for o(X) which is independent of 

the system’s parameters, calibration constants, C(X), are obtained with no 

gas in absorption cell by taking the ratio of the "through cell" signal to 

the "bypass cell" signal. When T (t) * 1 * T (f) and the dc 0 £f set * 
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The absorption cell is constructed with a cold finger as shown in 
figure 7. During the calibration run to determine C(.\), the cold finger 
is immersed in liquid nitrogen which condenses all the S0 2 in the cell. 
Following the calibration run, the liquid nitrogen is removed and the cold 
finger is allowed to return to room temperature. This allows the S0 2 to 
return to the path of the UV laser beam which will make T < 1. 

To determine the transmittance of the gas within the cell, the 
constants determined by equations (10) and (11) are multiplied by the ratio 
of the bypass cell signal to the through cell signal. Expressed 
mathematically. 


T (t , X) 

g 


c 0 ff U»X) 


Vu(t) 

vTrtT 


cell thawed (12) 


and 
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cell thawed (15) 


y*'*” 

The absorption coefficients for the off-line and on-line are then determined 
using equation (9) , and the wavelength is determined from a plot or table 
of absorption coefficients versus wavelength. The above analysis is based 
upon the assumptions that (1) the ^ c 0 ff set * 0 and remains 0 during the 
experiment, and (2) the wavelength at the time of determining C(X) is the 
same as when measuring T^(A), 


3,4. Timing and Control 

The timing and control are illustrated in figure 8 along with the out- 
puts from the photodiodes installed within the Nd-Yag laser and S0 2 absorp- 
tion cell assembly. The firing rate is controlled by a Berkeley Nucleonics 
Corp (BNC) model 8010 pulse generator #1. The complement positive output 
pulse is supplied to the EXTERNAL input jack of the ILS HL-101 laser and 
commands the laser to fire on its negative transition (r x ) . Logic and 
delays within the ILS laser result in the following sequence of signals 
with At set to 90 us and the DIP switches, or the flashlamp to Pockels 
cell command delay set to binary 61, which is 97 us: 


ti + 2 us 
t\ * 92 us 
tj + 100 us 

tj + 100 us 
tj + 102 us 
tj + 190 us 

* 192 us 


command to fire all three flashlamps 

command to fire all three flashlamps a second time 

leading edge (positive transition) of 50 us pulse 
produced on BNC output jack labeled MONITOR PULSE 

1 » t - 2 us 

o 

command to open Pockels cell and laser 
laser pulse A transmitted * t 

leading edge (positive transition) of 30 us pulse 
produced on BNC output jack labeled MONITOR PULSE 

2 * t ♦ 88 us 

o 

laser pulse B transmitted = t 0 + At * t + 90 us 


The monitor pulse #1 is used to trigger pulse generator #3 which outputs a 
]00-ns wide gate pulse with a transition from -0.5 V to -1.5 V required to 
be compatible with the ECl (Emitter Coupled Logic) used within the LeCroy 


r 


T 


Research Systems (LRS) model 227 quad integrator unit. This negative pulse 
is used to "gate on" channels i and 2 of the quad integrator starting at 
(t * 40 ns) and ending at (t Q + 60 ns). This will allow for a maximum of 
±40-ns jitter* of the 20-ns wide t pulse within the 100-ns gate. One 
hundred ns is the maximum gate duration recommended by the quad integrator 
manufacturer in order to obtain the specified linearity (±0,4% of reading, 
±0.2% of full scale) of integration. In addition, a narrow gate minimizes 
errors in the output due to the integration of noise on the signal line 
before and after the t pulse. 

The monitor pulse #2 triggers pulse generator t 4, which produces a 
100-ns wide gate pulse used to gate on channels 3 and 4 of the quad inte- 
grator. The delay time of this gate is set by front panel controls on 
pulse generator #4 starting at (t ♦ At - 40 ns) and ending at (t Q ♦ At + 

60 ns). 

The quad integrator outputs four voltages to the computer: Vi and V 2 

are the integrated t outputs from the through and pass absorption cell 
photodiodes respectively, and V3 and V** are the (t ♦ At) outputs from the 
through and pass absorption cell photodiodes respectively. These signals 
are processed by the computer as discussed in "Absorption Cell Theory." 

Pulse generator #2 is provided to furnish a Reset signal to the quad 
integrator and provide an advance pulse to the computer through a line 
driver. Both of these early pulses are provided at time t\. In addition, 
the computer is provided t and (t Q + At) through a signal conditioning 
and line driver module. These pulses are generated by a photodiode sensing 
the 532-nm laser pwwer envelope installed within the ILS laser pump. 

A voltage converter is coupled to pulse generator #1 to provide an 
18-V pulse to trigger a Holobeam model H-ll Pockels Cell Q-switch control 
unit. This unit produces a -13 ICV pulse starting at time t 2 . which is 
adjusted (by a potentiometer on the control unit) so that the Pockels cell 
(part of the optical switch) is fully open at time (t Q + At) when switch 
S-l is in the B on-line position. When switch S-l, located on the Pockets 
cell control unit, is placed in the A on-line position, the -13 KV pulse 


* The jitter was measured to be less than ±30 ns. 
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is repositioned ahead of t so that the Pockels cell is fully open at time 

t . Switch S-l then directs pulse A to either the off-line or the on-line 
o 

dye laser. 

A block diagram of the receiver system is shown in figure 9. The 
receiver consists of a 30-cm Cassegrain telescope with UV-enhanced coatings 
on the primary and secondary mirrors, and a detector package. In the 
detector package a dichroic beamsplitter is used to separate the back- 
scattered UV radiation near 300 nm from the backscattered visible radiation 
at 532 nm or 600 nm. The UV radiation is detected using an RCA C7268 photo 
multiplier tube which has the same dynode structure as an RCA 7265 photo- 
multiplier but with a bialkili photocathode for enhanced UV quantum effi- 
ciency. The visible radiation is detected using a standard RCA 7265 photo- 
multiplier tube. The gain of the photomultiplier tubes can be modulated so 
that it is proportional to t 2 , where t is time (ref. 11). This compensates 
for the dependence of the return signal and reduces the dynamic range 
of the signal. 

The signals from the photomultiplier tubes are fed to Biomation wave 
form recorders where they are converted to digital form. The digitised 
return signals are then read by the PDP 11/10 minicomputer and recorded on 
magnetic tape for later processing. The PDP 11/10 is programmed so that 
data analysis can proceed simultaneously with data recording. This permits 
preliminary analysis and display of the data in nearly real time so that 
operation of the system can be monitored and changed if necessary. The 
data acquisition system is the same as that used for the plume dispersion 
study and is described in more detail in Section 5. 

4. IR DIAL SYSTEM 

The infrared region of the spectrum from 1.5 to 4.5 vim contains funda- 
mental, combination, and overtone bands of a number of tropospheric pollu- 
tant and trace gases. The objective of this program is to develop a mobile 
DIAL system which would have the capability of making remote measurements 
of many of these gases with particular emphasis on effluents in power plant 
plumes. 
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Development of the system has been delayed by difficulties encountered 
by the laser supplier. However , progress has been made in assembling and 
testing other components of the system and in analytical studies to deter- 
mine operating wavelengths and expected sensitivities for some of the gases 
of interest. 


The most critical component of the system is the laser transmitter 
because it must have high power, tunability, and narrow line width. The 
laser is being developed by Westinghouse Research Laboratories under 
contract to NASA Langley Research Center, and is based on technology demon- 
strated at Stanford University (ref, 12) . A schematic of the laser trans- 
mitter is shown in figure 10. A Nd:YAG oscillator and amplifier is used to 
pump an Optical parametric Oscillator (OPO) and Optical Parametric Amplifier 
(OPA). The oscillator produces 150 mJ/pulse at 10 Hz with 0.1-cm" 1 line 
width from an unstable resonator cavity. The amplifier then increases the 
oscillator energy to 350 raJ/pulse. Approximately 60 mJ of the 1.06 pm 
energy are used to pump the OPO, and the remaining energy is used to pump 
the OPA. Angle tuned LiNb0 3 crystals are used in the OPO and OPA for the 
nonlinear conversion of 1.06-um energy to tunable IR energy. One of the 
requirements for wavelength conversion is + uk, where is the 

pump frequency, and u) and uk are the signal and idler frequencies, 
respectively. The signal wavelength is oscillated in the OPO cavity, and 
line narrowing is accomplished by a grating and etalon to achieve a line 
width in both signal and idler wavelengths of 0.1 cm* 1 . The OPA increases 

the output of the OPA at 2.2 ym to a total of 50 mJ/pulse in the signal and 

idler wavelengths. For DIAL measurements the OPO grating is rocked between 

modes of the etalon to produce alternating on/off signal and idler wave- 

lengths. Discrimination between the signal and idler wavelengths is accom- 
plished by filters in the receiver. 


The receiver system consists of a 46-cm diameter, gold-coated Newtonian 
telescope and liquid-nitrogen cooled, 1-mm diameter InSb detector. The 
laser beam must be transmitted coaxially with the telescope since the laser 
divergence is less than 0.4 mrad and the telescope field of view as deter- 
mined by the detector size is 0.8 mrad. This is accomplished using a 
Coude mount for the telescope so that the laser remains stationary and the 
laser beam is directed along the axes of rotation of the mount and then 
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coaxial with the telescope. An absorption cell with path length adjustable 
up to 200 m will be used to measure the absorption coefficient of the gas 
being measured for each laser pulse. Detectors for use with the absorption 
cell have been tested to determine their linear region of operation, and a 
gas-handling system for the absorption cell has been designed and built. 

The data acquisition system is similar to that used for the UV DIAL system. 

Available spectroscopic data have been used to make preliminary esti- 
mates of the measurement capability of the system for various gases. These 
estimates are shown in table 2 assuming a five percent uncertainty in meas- 
uring the transmittance. Improved sensitivity could be obtained by averaging 
laser shots. The system should have fairly high sensitivity for measuring 
HC1, CH 4 and H 2 O, The relative insensitivity of the system for measuring 
S0 2 limits its application to near stack exit conditions. The UV DIAL 
system is restricted to measuring concentrations of S0 2 which are less than 
300 ppm over a 10-m distance because of residual absorption by S0 2 at the 
off-wavelength, and therefore will be applied to measurements of S0 2 at 
large distances from stacks and in regional investigations. The lower 
sensitivity of the IR DIAL system thus complements the UV DIAL system in the 
measurement of S0 2 . 


Table 2. Estimated IR DIAL Sensitivities, 


Species 

Wavelength 

wm 

Optical Depth Sensitivity 
ppm - meter 

S0 2 

4.0 

540 

CO 

2.33 

660 

HC1 

3.S4 

25 

CHt, 

3.27 

5.3 

co 2 

2.06 

7400 

H 2 0 

Selectable 

Selectable 

n 2 o 

3.88 

280 
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5. POWER PLANT PLUME DISPERSION STUDY 


In early 1976 an agreement was made by NASA/LaRC and the Power Plant 
Siting Program of the State of Maryland Department of Natural Resources 
(PPSP) to perform a joint experiment to determine the utility of the lidar 
as a method for characterizing plume rise and plume dispersion from large 
fossil fuel burning power plants. The experiment was conducted during four 
one* week periods during September and October 1976. In this section the 
design and results of that experiment will be described. Some of the 
material in this section has been presented before (refs. 1, 4); however, 
it will be repeated here for completeness. 

The Power Plant Siting Program is responsible for determining the 
environmental impact of proposed power plants and for maintaining an inven* 
tory of environmentally suitable sites which can be acquired by a utility 
or exchanged for a site a utility owns which has been determined to be 
environmentally unacceptable. As an aid to improving the validity of 
environmental impact statements, PPSP has contracted with Martin Marietta 
Research Laboratories (ML) and Environmental Measurements Inc. (EMI) to 
perform a fairly extensive theoretical and experimental study of plume 
dispersion from tall stacks with the objective of improving existing plume 
rise and plume dispersion models, or developing new models if existing 
models prove to be inadequate. 

There were a number of organisations in addition to NASA/LaRC and 
PPSP who made contributions to the joint experiment. The experiment was 
conducted at the Morgantown Generating Plant of the Potomac Electric Power 
Company. Personnel at the plant were very cooperative in terms of providing 
information on operating conditions of the plant and providing space and 
facilities for making some of the measurements. MML served as a contact 

4 

with the power company and PPSP and coordinated the activities of the various 
groups in the field, EMI performed a major role in the experiment, making 
meteorological measurements and conventional plume measurements. Radio- 
sondes were launched three times a day to measure the temperature profile. 

In addition, pilot balloons were released every hour to measure wind speed 
and direction as a function of altitude, EMI plume measurements were made 
using two instrumented vans: one, the Air Quality Moving Laboratory (AQML) , 
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measured total integrated SO 2 overhead using a Barringer correlation spec- 
trometer and in situ SO 2 concentration as it traversed beneath the plume on 
available roads. The second van, the Transportable Laboratory (TPTL) was 
positioned near the plume centerline in an area where the plume was touching 
down and measured in situ SO 2 concentration on a time-averaged basis. 
Engineering-Test Services (ETS) made in- stack measurements of the mass emis- 
sion rate, SO 2 , total particulates, and particle size distribution. Lidar 
measurements of plume rise, plume dispersion, and mixing layer height were 
made by Old Dominion University and Wyle Laboratories. Lidar data reduction 
and analysis were performed by ODU using the NASA/LaRC computer facilities. 

The joint experiment was desired by a group which included the inves- 
tigators making the measurements and modelers and data analysts from NASA/ 
LaRC, MML, Goddard Institute for Space Studies, and ODU. The results of a 
simulation of the lidar experiment performed by Kibler (ref. 13) were used 
as a guide for selecting suitable sites for lidar operation and for devel- 
oping the optimum lidar data taking strategy. 

5.1. Mobile Lidar System 

The mobile lidar system was built by personnel from ODU, NASA, and 
Wyle Laboratories using NASA equipment. A block diagram of the system is 
shown in figure 11. The lidar system consists of a ruby laser, a tele- 
scope receiver, a detector package, and associated instrumentation. The 
laser, telescope, and detector package are mounted on a searchlight mount 
with tracking capability. The entire system, including instrumentation and 
the searchlight mount, is contained on a flatbed trailer. 

The laser used in the system is a Ho lobe am 600 Q- switched ruby laser 
with a beam divergence of 3*j mrad. The beam divergence was reduced to 
approximately 1 mrad by using an up-collimating telescope at the output of 
the laser. The output energy of the laser can vary from 0.75 to 2.0 J per 
pulse with a pulse length of 30 ns and a pulse repetition frequency up to 
1 pulse per second. The ultimate range resolution attainable is 4.5 m. 

The receiver is a 30-cm cassegrain type telescope with a field of 
view of approximately 4 mrad. The detector package has a provision for 
mounting two photomultiplier tubes for extended dynamic range. For the 
plume dispersion experiment, only one tube was used. An RCA 7265 photo- 
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multiplier tube was selected for use In the system since it is sensitive to 
the 694-nm laser radiation, and it could be easily gated to prevent overload 
from the close-in return. 

The instrumentation used consists of a high voltage power supply and 
gating circuit for the photomultiplier tube, pulse generators used for 
timing, and the data acquisition system. 

The data acquisition system is based on a Digital Equipment Corporation 
model POP 11/10 minicomputer. The lidar return signal from the photomulti- 
plier passes through a log amplifier and is recorded by a Biomat ion model 
8100 transient digitizer which can record 2048 8-bit words of data at sample 
rates up to 100 MHz. The log amplifier reduces the dynamic range of the 
signal to the transient digitizer. The transient digitizer has an analog 
output so that the data can be displayed continuously on an oscilloscope 
and a digital output from which the computer accepts the recorded data. 

Other data accepted by the computer include laser energy, laser shot 
counter reading, and elevation and azimuth angle from the searchlight mount. 
The data from each laser firing is immediately recorded on magnetic tape 
for later processing using either the PDP 11/10 or the main Langley Research 
Center computer facility. Also attached to the minicomputer is a Ramtek 
graphics display generator which can display information in 16 shades of 
gray on a standard black and white television monitor. 

The computer has been programed to allow data processing to proceed 
simultaneously with data recording. This permits preliminary data analysis 
(such as range correction) and display on a nearly real-time basis. Three 
different types of data display available: one is called an A-scope 

display and is simply an X-Y ith no intensity modulation of the 

display. Another type of disp.-, , called a 2-scope display, is useful for 
applications where the lidar is pointed vertically, such as in mixing layer 
height measurements. This display plots intensity versus height on a 
vertical line using the 16-shade gray scale modulation. The third type of 
display is the RHI or range, height, intensity display which is used for 
displaying plume dispersion lidar returns. For this display, the horizontal 
axis corresponds to the horizontal distance from the lidar and the vertical 
axis corresponds to height above the lidar. Intensity is then plotted 
along a line which corresponds to the lidar elevation. 
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The maximum range o *2 the lidar is not known. However, in this experi- 
ment measurements of plumes were made up to nine km from the power plant 
and up to nine km from the lidar. 

The Morgantown power plant is located on the Potomac River in Charles 
County, Maryland about 56 las (35 mi) south of Washington, D.C. where U. S. 
Route 301 crosses the Potomac River. The plant has two 575-MW generating 
units operating on coal or oil or a mixtuie of both. During the time when 
this experiment was conducted, the fuel mix was approximately half coal and 
half oil. Flue gases are exhausted through two 213-a stacks located 76 m 
apart. The stack exit velocity is 30-m/s and the stack exit temperature is 
13$*C, For the two one-week measurement periods in September, both gener- 
ating units were operating. However, for the two one-week periods in 
October, one unit was shut down for maintenance, so only one unit was 
operating. 

Figure 12 shows a map of the area around the power plant with the lidar 
sites used indicated by stars. Circles centered at the power plant with 
one-km increments of radius are drawn on the map for reference. 

The basic measurement strategy was to locate the lidar two to four km 
from the plant at a site selected so that the plume could be viewed from 
the side. Cross sections of the plume were obtained by scanning the lidar 
vertically through the plume at selected downwind distances from the stack. 
The measurement sequence was repeated several times over a one-hour period 
in order to obtain time-averaged measurements. Where possible, one of the 
downwind distances was selected to coincide with simultaneous measurements 
made by EMI using the AQML. Lidar measurements of mixing layer height were 
made before and 'fter each series of plume measurements. 

Because of the time required to relocate the lidar (approximately six 
hours) and the potential electronic malfunctions which could be caused by 
moving, a lidar site was chosen for several days of operation even though 
it might not be ideal for every day. The site chosen was based on the 
weather forecast from the National Weather Service. As it turned out. one 
site was used for each one-week measurement period. For the period from 
September 7-10, 197o, the lidar was located at site number 1 indicated in 
figure 12 at a distance of 1.91 km and bearing 74 s from the north stack of 
the power plant. From September 20-24 the lidar was located at site 


number 2, 3.12 km from the north stack at a bearing of 145 s . From October 
11-1S the lidar was located at site number 3, 3.84 km from the north stack 
at a bearing of 5.8*. From October 25*29 the lidar was located at site 
number 4, 3.54 ka from the north stack at a bearing of 261*. 

Data reduction and analysis were performed using the LaPC computer 
fr 'ility. 


5.2. Data Reduction Techniques of Lidar Plume Measurements 

The returns are initially corrected for a system offset and log aapli* 
fication. The offset is required due to inherent limitations of the log 
amplifier. Actual amplification differs from the true log function, so 
calibration techniques were used to determine a log transfer function. 

These corrections to the lidar signal are combined in the following 
equation: 

V* (R) « 5.8063 + 0.0337 V(R) - 0.00004 V 2 (R) (14) 

A background return is obtained by gating on the photomultiplier tube 
without firing the laser. This signal represents ambient light entering the 
tube as well as electronic noise. Since the background return is fed 
through the log amplifier it must al«o be corrected by equation (14). 

Once the plume and background returns are corrected, the bacKground 
is subtracted and the result is range corrected by R 2 where range (in km) is 
determined as follows: 

R. ■ 0.15 (t ♦ 1/f) i - 1,2024 

t ■ initial delay (usee) 
f ^ sampling frequency (MHz) (15) 

The range resolution for a sampling frequency of 50 Mir is then 5 ra. Alti- 
tude (2) and horizontal distance from the lidar (Y) are calculated using 
the elevatior. angle ( 0 ): 



Z » sin (6) + Zsi 
Y * cos (8) 

Zsl * altitude above sea level of lidar. 


( 16 ) 

(17) 


After correcting the raw data for log amplification, background and 
range, the resulting profiles represent relative measurements of the atten- 
uated aerosol and molecular atmosphere along the lidar line of sight. To 
eliminate the attenuation effects, a profile of the molecular atmosphere is 
usually employed to normalize the lidar returns. Since a molecular profile 
was not available, a portion of the lidar return, prior to the plum and 
extrapolated on out, was chosen to represent the ambient atmosphere. Due 
to spurious tails introduced by misalignment of the lidar system, a quadratic 
fit is used. The limits of the curve fitting are specified and are constant 
for any set of scans at one particular azimuth angle. The lidar profiles 
are then divided by the fitted curve to yield ratios of plume to nonplume 
data. The curve-fitting limits are chosen (by trial and error) to produce 
ratios close to unity for nonplume portions of the profile. In some cases, 
the curve had to be fitted to sections of the profile both before and after 
the plume to obtain these results. Truncation is necessary in some cases 
at one or both ends of the profiles. To obtain relative magnitudes of only 
the plume, one is subtracted from a. I ratios, and any results less than 
zero are set equal to zero. 

Grid limits are chosen so as to center the plume in both the Y and Z 
planes. Profiles are placed in the grid, and where several data points 
fall within the same grid box an average value is computed. Interpola- 
tion is performed on any empty grid box adjacent to a nonempty grid box 
using the four adjacent boxes. A negative sign is attached to interpolated 
points to differentiate them from actual values. The value of -999.0 is 
placed in all remaining empty boxes to indicate lack rf data rather than 
lack of plume. 

The total burden of a scan is calculated by summing all actual and 
interpolated points in the grid. If the lidar measure of the plum*' is 
denoted by S then the total burden (S^) is: 


(IS) 




s 
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The centroids (Y and I s ) and their deviations (o v and cj„) are defined by 

T «# 


(t ? s v >' s t 

(19) 

(? ? s ’“> s t 

(20) 

($ £ (W*}/S t - 

(21) 

(£ £ (S 1 *)/$^ - ?)'/5 

(22) 


The rotation correction (a) previously mentioned and the centroid distance 
downwind (X) are determined as follows: 


)C • ♦ Y* • L V cos ~ 1 (23) 

a » cos |*/2 * c°s^|^ * T 2 - 1^/(2 X Y)J^ (24) 

where L is the distance of lidar from the smokestack and 4 s is the azimuth 
angle of scan. 

Two types of averages were calculated for a set of scans at a given 
azimuth angle: an Fulerlan average and a Lagrangian average. For the 
Eulerlan plume average, grid points, including interpolated values, are 
averaged over a set of scans at the same azimuth angle. No further inter- 
polation is performed. The values of S . Y, Z, <?y» o„, Y and <\ are then 
recalculated for the averaged scan exact l v as before. 

For the Lagrangian plume average the individual grids are averaged 
together by superimposing the centroids. The new grid is related to the 
original grids by the following: 


Y* * (Y * T) a 

(25) 

:* * (Z - Z1 

(26) 



The values of interest are again recalculated in a similar manner as 
before with the following exceptions: 

(1) The rotation correction is eliminated from equation (19) 
since it has already been accounted for in equation (25). 

(2) The value of Y in equation (23) is an average of the Y 
centroids on the individual scans. The Y centroid on the 
Lagrangian average should be zero, but due to the shifting 
of the new grid this is not usually the case. 

(3) The rotation correction is not calculated in equation (24) 
since individual grids are rotated prior to placement in 
this averaged grid [eq. (25)]. 

Once the grids representing individual scans and averages of scans 
are created, a number of other calculations can be made and the results 
plotted or presented in tabular form. 

The grid is divided into 10 equally spaced vertical and horizontal 
slices. The following computations are then made: 

(1) Maximum S within the slice (plots 1 and 6) and its corresponding 
Y value (plot 2) or corresponding Z value (plot 7) . 

(2) Line integrals determined by summing all values of S within 
a horizontal slice (plot 3) or within a vertical slice 
(plot 8). 

(3) The variation of the centroid and its deviation versus altitude 
(plot 4) and versus distance from the lidar (plot 9). 

Plots 5 and 10 represent slices up the centroids of the plume and the 
normal curve computed from the mean and standard deviation. The slices on 
these two plots are of the same resolution as the grid, which at present 
is 20 ra. 

In addition, a contour plot is available for each individual and each 
averaged grid. 

A magnetic tape has been generated which contains all the lidar plume 
data obtained in the experiment. The data is recorded in the form of the 
grids described previously. 
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Each grid is contained in a separate file which consists of a header 
record, the grid, and end of file mark. The information in the header 
record is tabulated in table 3. The number of points in the grid are 
stored in words 22 and 25. The grid was placed on tape with the 
following implied do-loop: 

JTOT - HDR(22) (2 - plane) 

ITOT - HDR(2S) (Y - plane) 

((S(I,J), I « 1, ITOT), J ■ 1, JTOT) 

The entire tape is written with a 10E10.3 format. There are a total 
of 423 grids (either individual or averaged) on the tape from the Morgan- 
town study. 

A copy of the data tape and the plots and tables described, as well as 
the meteorological data and conventional plume measurements made by EMI, 
may be obtained by writing to: 

Dr. F. S. Mills 
Mail Stop 401A 

NASA Langley Research Center 
Hampton, Virginia 23665. 



Table 3. Header record for plume data (format 10E10.3), 


Word No. 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 


11 

12 

13 

14 

15 


month ) 
day > 
year ; 

hour 1 
minute i 
second ) 

hour ) 
minute V 
second ) 


date 


start time of scan 


end time of scan 


lidar site 1 ■ PEPCO 

2 * Waver ly Pt. 

3 * Watson's Farm 

4 * Dahlgren Beach 

total burden 

minimum elevation of scan 
maximum elevation of scan 
aximuth angle of scan 
sea level altitude at lidar site 


16 angle 

17 distance 


relative to north stack 


18 rotation correction 


19 

minimum Z 



20 

maximum Z 



21 

grid resolution 



22 

number of points 

in 

Z direction 

23 

minimum Y 



24 

maximum Y 



25 

number of points 

in 

Y direction 


grid dimensions 


26 file number where scan is placed on tape 

27 number of scans in grid 


28 Y centroid 

29 Z centroid 

30 centroid distance from stack 

31 Y sigma 

32 Z sigma 



ALTITUOE (KM) 



Figure 1. Calculated water vapor concentration using electronically 
simulated return signals. 
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Figure 2. Effects of two types of error on calculated 
concentration. 
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Figure 4. Double-pulse tunable UV dye laser. 
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Figure 8. UV DIAL lidar tiding and control using aonitors 1 and 2 froa ILS. 
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Figure 10. Block diagram of the IR tunable laser. 
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Figure 11. Block diagram of the plume dispersion Hdar system. 
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